Investigation of resonant and transient phenomena in Josephson junction flux qubits 
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We present an analytical and computational study of resonances and transient responses in a 
classical Josephson junction system. A theoretical basis for resonances in a superconducting loop 
with three junctions is presented, outlining both the direct relationship between the dynamics of 
single- and mult i-j unction systems, and the direct relationships between observations of the classical 
counterparts to Rabi oscillations, Ramsey fringes, and spin echo oscillations in this class of systems. 
We show simulation data along with analytical analyses of the classical model, and the results 
are related to previously reported experiments conducted on three junction loops. We further 
investigate the effect of off-resonant microwave perturbations to, e.g., the Rabi- type response of the 
Josephson system, and we relate this response back to the nonlinear and multi- valued resonance 
behavior previously reported for a single Josephson junction. The close relationships between single 
and multi-junction behavior demonstrates the underlying dynamical mechanism for a whole class of 
classical counterparts to expected quantum mechanical observations in a variety of systems; namely 
the resonant and transient behavior of a particle in an anharmonic potential well with subsequent 
escape. 

PACS numbers: 74.50.-l-r,03.67.Lx,85.25.Cp 

I. INTRODUCTION 

Since the mid 1980's several experiments have sought to establish that Josephson systems show evidence of quantum 
interference of macroscopically distinct states. Various configurations have been utilized to elicit the phenomenon, 
including current-biased single-junction circuits [1-5], flux-biased two-junction SQUIDs [6], charge-biased two-junction 
SQUIDs [7], flux-biased three-junction SQUIDs [8-13], charge-biased three-junction SQUIDs [14], and hybrid systems 
which often include inductor-based systems to create multiple potential wells [15-19]. The phenomena observed to 
date include multi-peaked probability distributions in bias current, Rabi oscillations [20], Ramsey fringes [21], spin-echo 
oscillations [22], and quantum tomography[23]. 

In several recent investigations the authors of the present paper have contributed to an alternative interpretation of 
the basic types of experiments on Josephson junctions. This interpretation relies on the well established phenomenol- 
ogy of the Resistively and Capacitively Shunted Junction (RCSJ) model [24] and on nonlinear dynamics related to 
it. Multi-peaked resonances were addressed in Refs. [25-27] while References [28,29] developed analytical expressions 
and presented computer simulations which described the relationships between classical Rabi-type oscillations, driv- 
ing amplitude and frequency, and systematic dissipation. The related phenomena of Ramsey-fringes and spin-echo 
oscillations are depicted and described from a classical point of view in Refs. [30,31]. While this exemplifies the 
phenomena through current-biased single-junction circuits, we submit that a wide class of Josephson systems exhibits 
same behavior when a system variable is trapped in a potential well [32]. This presentation will amplify this point. 

A few examples of multi-junction systems for which direct comparisons between classical theory and experimental 
results have been made in the low temperature regime, where quantum phenomena should be observed, include multi- 
peaked resonances in Josephson interferometers, made of two junctions connected in parallel [33], and a more complex 
Josephson interferometer-based device[34] containing a double-SQUID, in which a butterfly catastrophe was observed 
at 10 mK. In both cases, the standard classical RCSJ model provided sufficient detail to understand the experimental 
observations. 

The observations from interferometer systems described in the last paragraph encouraged us [32] to extend the 
classical ideas and implementations developed for the single junctions and two-junction interferometers to the more 
complex system of the fiux-biased three-junction loop [10-12, 35]. It is the plan of the present paper to report on our 
approach to this problem. In this approach the analysis of the most important features of three junctions loop under 
the influence of pulsed microwaves was preluded by a study of analogous single junction phenomenology. 

In the next section we derive the basic equations for the three junction loop and present an analytical description of 
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resonances in the system. In sect. 3 wc demonstrate computer simulations which map the resonances as a function of 
the biases (the magnetic flux in the loop). Simulations are shown which depict the Rabi-type oscillations, Ramsey- type 
fringes and spin-echo-type oscillations in the three-junction loop. Also, off-resonant behavior in the three-junction 
loop is explored. 

Our simulations and analytical treatment allow us to interpret the essential experimental features reported for 
Josephson circuits; the analysis shall all be based on the RCSJ modelling of Josephon circuits[24] who took origin 
from early remarks on the classical (i.e. , continuous) nature of Josephson phase and flux variables. In this context 
we also refer to the subjects of our presentation as "classical" analyses and modelling. 



II. THREE-JUNCTION LOOP RCSJ MODEL 



A. Single-junction theory 

Figure 1(a) depicts the single Josephson junction RCSJ model. Such a physical system is also referred to as 
Josephson phase qubit. The normalized classical equation can be writtcn[24] 

(p + aip + sinip = r] + Ss sm{u>gt + Os) + ep{t) + n{t), (1) 

where ip is the difference between the phases of the quantum mechanical wave functions defining the junction, rj 
represents the dc bias current, and es(t) and LOg represent microwave ciirrcnt amplitude and frequency, respectively. 
All currents are normalized to the critical current /c, and time is measured in units of the inverse plasma frequency, 
LVo~^, where wq^ = 2eIc/hC = 2nIc/^oC, C being the capacitance of the junction and $o = h/2e the flux quantum. 
Tunneling of quasiparticles is represented by the dissipative term, where a = HuJo/2eRIc is given by the shunt 
resistance R, and the accompanying thermal fluctuations are time-correlated by the normalized fluctuation-dissipation 
relationship [36] 

(n(i)) = (2) 
(n(i)n(i')) = 2a'^6{t-t')=2aeS{t-t'), (3) 

T being the temperature and Hj is the characteristic Josephson energy Hj — Ich/2e. 
The normalized energy, {H — Hq) is deflned by 

H = + 1 — cos (p — rjif (4) 

Ho = 1 - a/I - - r^sin^^ 77 , (5) 

where Hq is the energy at the minimum of a potential well. In Refs. [28,29] a classical analysis is presented for 
the modulated frequency, which is small when compared to the driving signal ujs, and referred to as the Rabi-type 
oscillation frequency Qji due to its direct relationship with experimentally reported Rabi oscillations. Thus, the 
classical analog to quantum mechanical Rabi oscillations has been denoted Rabi-type oscillations [28]. 



B. Three junction loop theory 



Figure lb depicts the three-junction loop which will be the subject of the analysis of this paper. We concentrate 
our study on this system which has been the core of several experiments and theoretical analyses [35]. It consists of 
a superconducting loop containing two identical junctions and one with smaller capacity and critical current. The 
normalized classical equations can be written [35, 37], 

(pi + aipi + sinipi = -J-[2nM{t) + (pi + ip2 + ipa] + ni{t) 

PL 

ip2 + aip2 + sin ip2 = — [2TrM{t) + ipi + ip2 + Va] + n2{t) (6) 

PL 

03 + Oi(f>z -\- sin ips = —5- [2iTM{t) + ip>i+ip2 + Vs] + ns{t). 
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FIG. 1: Circuit diagram for (a) the single junction model and (b) the three junction loop. In (b) the smaller junction is 
indicated by the factor k appklied to the capacitance, resistance , and critical current. . 



where ifi is the difference between the phases of the quantum mechanical wave functions defining the ith junction, k 
is the scale factor relating the one smaller junction, represented by (^3, to the other two larger junctions, /3l is the 
normalized loop inductance /J^ = 2nLIc/^o. M{t) is a variable which summes up the external magnetic fluxes to 
which the system can be exposed 

M{t) = Mdc + es{t) sin(w,t + 6*,) + ep{t), (7) 

where Mdc represents the flux of a time-independent fleld, £s{t) and tOg represent the ac flux (signal) amplitude and 
frequency, respectively while a flux pulse for probing the state of the system is here taken into account by £p{t). A 
general analysis of the above system of equations for arbitrary values of the parameters would surely lead to rather 
intriguing dynamics and complex oscillations. In the analysis that we will present herein we will focus on the limits 
that are suggested by the experimental reality. A very indicative example of the effects that a flnite inductance can 
have on Josephson flux qubit was reported in ref. 34. 

We write the Hamiltonian for the undamped {a = 0) system as 

H = Hk + Hp 

Hp — 2 ~ cos Lpi — cos ip2 + k{1 

For small Pl in Eq. (8), we can see that energy from the loop current (last term) dominates the sum of the phases. 
We introduce now a transformation which allows combining the phases ipi, </32,and 933 while maintaining the kinetic 
part of the Hamiltonian still separable in three distinct variables (jPit^2, and "03) 



(8) 



cos (^3) + 77-^('^i + (P2+ 'P3 + 27rM)^ 
2pL 



(9) 
(10) 



tpl = ipi+(p2 + V3 ^ ( fl = TWV'I + 5^2 - ^03 

02 = - '/'2 > I f2 = J^lpl - \i>2 - |03 

^3 = -TT2l(¥'l + ^2) + J^^3 J [^3 = YI^V-l + 03 
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The Hamiltonian now becomes 
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"01 + V'3 



-^(Vi+27rM)2. 



Under the condition <^ 1, this system can be simphfied to a single degree of freedom by examining the constraints 
which frame V^i and ^2- First, it can be seen from Eq. (11) that potential energy considerations prevent the deviation 
of ■^i far from —2-kM. 

A second constraint, on ^2, is obtained by considering the equation of motion 



V'2 + a'4}2 + 2 cos 



1 



1 + 2k' 



V'l - oV's sin — = ni(t) - n2{t). 



Prom this it can seen that the static equilibrium value of -02 is given by 



cos ■ 



i'2 



1 for cos(Yq^Vi - 5^3) > 
-I for cos(Yq^Vi - ^^'3) < 



(12) 



Note that this is a weak constraint as it does not require much energy to break the equilibrium. Equation (11) can 
now be written 



K 1 1 

Hp = 2 T 2 cos( ^^^J ^i- 2^3) + «(1 - cos(y 
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(13) 



where ± refers to the optimized choice based on Eq. (12). The constraint on ■^i can be determined by minimizing Hp 
with respect to ■01 



dHp 
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which is correct to first order in Pl- 

The potential energy can now be aproximated by a function of V'3 only 
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where the choice between =F should be such that the energy is minimized. 

Figure 2 depicts this single-degree-of-freedom representation of the potential energy. Note that while, in general, 
this potential consists of pairs of wells separated by higher energy "cusps" , for purposes of this work we focus on a 
single well-pair, as the energies are kept below these cusp 

Calculating the linear resonance frequency for a given well, the equation of motion for the friction-less and constant- 
M system is 



1, , •• dHr, 

- 1 + 2Kh^3 + ^ = 0, 

2 dips 



such that the fixed point V3 is given by 
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FIG. 2: Three-junction loop Potential Energy in terms of the reduced variable V'3- Each view represents a different value of 
Mdc: (a) Mdc = 0.5, (b) Mdc = 0.52, and (c) Mdc = 0.542. Other parameters are k = 0.68 and /3 = 0.09559. 



and for small oscillations 6 around the fixed point ips 



+ 6 provide the linear resonance frequency oji 
6 = 0^ 



I2 ^^^^ 



9^1 



/(1 + 2k) . 



(14) 
(15) 



In the next section we will use this result to characterize the switching statistics of microwave induced excitations in 
the classical model of the three-junction superconducting loop. 



III. SIMULATION RESULTS 



Our system is periodic under the application of the extrenal magnetic flux. In order to best describe the experimental 
reality the temperatures that wc consider arc in the tens of millikclvin range: we point out, however, that the kind of 
phenomena like those that we shall describe herein have been characterized in a temperature range spanning from few 
tens of mK up to 800 mK[32]; we believe that most of the phenomena herein described would display quite similar 
dependencies on temperature over the same range. 

To verify equation (14), we conducted computer simulations of pulsed microwave spectroscopy for the three-junction 
loop described by equations (6); naturally, the simulations are performed on the entire system of equations regardless 
of the analytical approximations. Microwave spectroscopy, combined with thermal escape measurements has turned 
out to be a reliable tool for investing resonances[25-27]; here we try to follow real experimental recipes[33]; we must 
bear in mind that for the thermal escape we simulate the whole system of Eq.6 and that switches to voltage states 
are just detected across a parallel connection of two junctions which forces us to attribute a 50% equal probability of 
switching to the two junctions. The results of our tracing the resonance through thermal escape and microvave pulses 
arc shown in Fig. 3. Each open circle represents the minimum flux pulse Sp required to obtain a statistical response 
of 50% escape probability. Each trial in the simulation is conducted in the following manner: with initial conditions 
of V'3 in the minimum energy configuration and Vs = 0, the system is driven at frequency uig and amplitude Sg for 
a period of {2 /a). At the end of the driving phase and after a slight pause (50 time units), a pulse of magnitude 
Ep is applied (in similar fashion to reference [30]) and the subsequent escape or non-escape to the alternate well is 
recorded. Statistics are then gathered to determine Sp^'^" for varying frequency. A resonance is determined for each 
value of Mdc, as indicated by a minimum value of Sp^^" as a function of frequency. In the inset of Fig. 3 we show the 
resonance effect in terms of switching probability, exemplified 

The solid line in Fig. 3 represents the resonant frequency response as predicted by Eq. (14). The open circles 
indicate the simulation results and depict the minimum amplitude probe pulse Sp for which the escape rate equals 
50%. The inset shows the relationship between probe pulse (at 50% escape probability) and signal frequency for 
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FIG. 3; Resonant frequency response in the three-junction loop. Lines represent the predicted frequency response from Eq. (14). 
The open circles indicate the resonances from the results of our sinmlations. The inset is a switching distribution indicating the 
resonance at the box marked 'A'. This provides an example of the placement for the open circles in the larger plot. Parameters 
for the simulation were Pl = 0.09994, k = 0.68, a = 4 • 10"®, = 4 • 10"*, and G(T) = 9.12 • 10"*. The dashed line indicates 
the symmetric relationship in linear resonance frequency. The data presented is based on 2500-7000 escape events. 



Mrfc = 0.48. The minimuin energy value in each well is marked with a boxed symbol to establish the correspondence 

with the potential energy plot. The dashed line indicates the symmetric relationship in linear resonance frequency 
(about Mdc — 0.5) for a well placement convention which is opposite to that of the continuous curve: in other terms 
the continuous line we trace the resonance moving the lower energy well (follow Fig. 2) from left to right while the 
dashed line would correspond to tracing the resonance moving the lower energy well from right to left. The data 
presented is based on 2500-7000 escape events ; the simulation results show close agreement with our theory although 
the agreement diminishes slightly for increasing values of Mdc- This can be understood by noting that as the flux bias 
is increased, the higher of the two wells becomes broader, allowing larger oscillations, which corresponds to a greater 
degree of anharmonic behavior. 

Next wc investigated whether the results obtained on the linear resonance freqiiency calciilatecl from Eq. 14 will 
generate features consistent with what one could expect from spectroscopic measurements on real Josephson junction 
systems. Thus we tuned the parameters of our simulations to generate Rabi-type oscillations, Ramsey-type fringes, and 
spin-echo- type oscillations [30, 31] [32]. The recipes for generating these oscillations follow closely the ones described in 
the previous publications for the single junction model and a summary of the signalling used for generating Rabi-type 
oscillations, Ramsey-type fringes, and spin-echo is shown in Figure 4 together with the idealized phase responses in 
the two relevant cases of escape and no-escape from the potential well. Our statistical simulation results are shown 
in Figs. 5 and 6 . The normalized slope that we extract from the linear dependency of the Rabi frequency versus ac 
amplitude of Fig. 5b is very close to that of the experimental results reported in Fig. 4b of Ref. [12] and with Fig. 
4b of Ref. [10]. We conclude that our numerical and analytical results are quite consistent with the experimental 
reality. As far as the Ramsey-type fringes are concerned, if we compare Fig. 5d in this presentation with the results 
reported in Fig. 4d of Ref. [12], we also find very good agreement in the Ramsey-type fringe being represented 
by a unity-sloped V-shaped fringe frequency dependence on the microwave detuning from the oscillator resonance. 
We clearly see our V-shape following the dashed lines in Fig. 5d. Thus, just as it has been demonstrated for Rabi 
oscillations, the classical system exhibits a signature similar, if not identical, to what one would expect from quantum 
theory. 

We note that our analytical approximations explain how the three-loop potential can be reduced to a single degree 
of freedom potential, however, it is clear from the results of the direct integration of eqs. (6) shown in Fig. 6 that the 
response of the loop is quite similar to that of a single junction [30, 31]. 




time ( arbitrary units ) 

FIG. 4: Sketches of idealized microwave perturbations and corresponding classical system (phase) responses for (a) Rabi 
oscillations, (b) Ramsey fringes, and (c) spin echo. Upper part of each plot illustrates a situation in which the probe (read out) 
pulse results in an escape from the potential well, while the lower part illustrates a no-escape case. The experimentally varied 
time for observing the oscillations in response is shown on each plot as At. 



The classical result, leading to the Ramsey fringe frequency being identical to the detuning between the applied 
microwave frequency u)s and the intrinsic resonance cJr of the oscillator, can be rationalized by considering the system 
behavior in the interval between the two 7r/2-pulscs. The oscillator is, at the time of the conclusion of the first 7r/2- 
pulse, defined by a specific phase-relationship (phase-locking) to the phase of the 7r-pulse. This phase-relationship is 
free to detune in the (ballistic) interval At between the two 7r/2-pulses, and the mutual phase detuning 6f between 
the microwave field and the oscillator is therefore given by 

/■At 

dp = I (Ws -u)r{t))dt , 
Jo 

where Wr can be a (weak) function of time due to the anharmonicity of the potential energy; i.e., as the oscillation 
amplitude decreases, due to the intrinsic damping in the system, the natural resonance frequency of the oscillator 
increases. However, for a lightly damped system, where aAt <C 1, the above time integration can be simplified, and 
the result for the Ramsey fringe frequency becomes 

o 

S iiT = — W UJg -U)r , 

which is exactly the detuning between the applied microwave frequency and the natural resonance frequency of the 

oscillator in the phase-locked state. Thus, the observed Ramsey fringe frequency dependence on the microwave 
detuning from resonance is easily understood classically for this system, and the intrinsic oscillator frequency ojr < u>i 
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FIG. 5: Oscillation phenomena for the three-junction loop: Rabi-type oscillations and Ramsey-type fringes. Panel (a) shows 
Rabi-type oscillations. Panel (b) indicates the Rabi-type frequency as a function of microwave amplitude, with Ss for lJs = 
0.45629 ~ LOr- Panel (c) is the resulting switching distribution for Ramsey-type fringes with ujs = 0.42528 and Sp = 0.0195. The 
inset depicts the driving frequency which achieves a fringe frequency of zero; here = 0.45629 ~ Wr, Ep = 0.0148. Panel (d) 
provides the relationship between fringe frequency and driving frequency. The arrow indicates measurement of uji = 0.46557 
by direct simulation. Note Eq. (14) predicts u>i = 0.46625. Parameters for all panels, unless otherwise noted, were M^c = 0.52, 
Pl = 0.09559, K = 0.68, a = 1.5 • 10"*, = 1.82 • 10"^, ujs = 0.443, Sp = 0.0149, and B{T) = 3 • 10"^; statistics were gathered 
for ~ 20 000 escape events. 



is close to the microwave frequency that results in fJ^^ 0. We submit that this detuning is also causing the so- 
called spin-echo measurements, which are generated by maintaining a constant temporal separation between the two 
7r/2-pulses from Ramsey fringe measurements, and then inserting a 7r-pulse. The final phase-relationship between the 
oscillator and the external microwave signal is detected by the second 7r/2-pulse, but this phase-relationship is now 
subject to the phase-twisting caused by the inserted 7r-pulse, which in turn depends on the detuning at the time of 
TT-pulse initiation. Thus, while the Rabi-frequency determines the magnitude and duration of the microwave pulses 
used in Ramsey- fringe experiments, it is the Ramsey- fringe frequency that determines the echo from the detuning at 
the time of the initiation of the 7r-pulse [31]. 

In previous work, and here in Fig. 5b, we have considered Rabi-type oscillations using the resonant frequency as 
the driving frequency: = ~ w^- However, this is not necessarily the case for all experiments. In Fig. 7 we show 
the effect of a lower-than-resonance {uig < coi) driving signal on Rabi-type frequency. One important aspect of these 
plots is that low-amplitude signals result in modulation frequencies which do not lie on the "main sequence" curve 
to which higher-amplitude responses conform. We assert that the reason for this is found in the multi-valued nature 
of the amplitude (and energy) response for off-resonant driving as shown in Ref. [27], Fig. 1. Also notable is the fact 
that this effect (an analytical result for the single junction) is also seen in the three-junction loop. 

IV. CONCLUSIONS 

The analysis we have developed reducing the three-junction loop to a single degree of freedom provides a significant 
reduction in the complexity of the system. The resulting equation provides direct agreement with our simulations 
when identifying the system resonances from experimentally relevant switching experiments. 

Our studies of off-resonant driving of Rabi-type oscillations reinforce earlier work concerning multi- valued functions 
in signal amplitude. Although no distinct conclusions can be made with regard to a connection between off-resonant 
driving signals and linear relationships between Rabi-type frequency and signal amplitude, we have shown evidence 
that driving the system near (but not exactly at) resonance does extract nearly linear behavior in the response for 
larger amplitudes. We emphasize that the simple model we are using may not completely account for all the details 
of the observed phenomena. 
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FIG. 6: Spin-echo-type oscillations for three-junction loop. Parameters are Mdc = 0.52, fiL = 0.09559, k = 0.68, a = 1.5- 10 
ujs = 0.419787, Es = 1.82 • 10"^, Sj, = 0.0195, and e(T) = 3 • 10"^ Each dot represents 22000 escape events. 
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FIG. 7: Simulations if off-resonance Rabi-type oscillation frequencies. Panel (a) depicts the off-resonance modes for the three- 
junction loop. Parameters were: Mdc = 0.52, Pl = 0.09559, k = 0.68, ujs = 0.443, = 1.82-10"^, Sp = 0.063, and 6 = 3-10"^. 
Two values of characteristic damping are given. The filled squares correspond to a = 1.5 • 10~*, empty circles:a = 0. The 
resonance frequency, u>r (as shown in Fig. 3d) is 0.45629. Panel (b) shows the off-resonance modes for the single-junction 
circuit. Parameters were a ~ 1.5 • 10"*, r? = yjl-u)f, ujs = 0.99lji = 0.646188, and 6 = 2 • 10"". 



Consistently with our previous work, wo have used the classical RCSJ phenomenology to investigate various phe- 
nomena heretofore attributed to macroscopic quantum tunneling and we have shown very good agreement with 
experimental results so far reported. We have shown that, with regard to Rabi-type oscillations, Ramsey-type fringes, 
and spin-echo-type oscillations, the three-junction loop differs little, qualitatively, from the single-junction circuit. 
This is easily attributed to the fact that the important dynamics of these systems is governed by resonant excitations 
in anharmonic potential wells followed by a read-out perturbation which makes the system escape from the well. Thus, 
our developed analysis and intuition about the system response to applied microwave pulses translates easily from a 
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simple single junction system to other systems with similar overall resonant excitations and escape from anharmonic 
wells. The main difference between the single junction system and other systems, such as the three loop system stud- 
ied here, is that the experimental reality may be different. However, given the relative insensitivity to noise inherent 
in the three-junction loop, it is understandable that several research groups work with flux qubits rather than single 
junctions (biased by dc bias current supplies). 

The work herein presented was developed on the background of the PhD thesis [32] of one of the authors (JEM) 
which was submitted to the University of California at Davis in the fall of 2006. In the present paper we have chosen 
to keep the same terminolology (Rabi-type oscillations, Ramsey-type fringes etc) introduced in the thesis and in the 
papers which inspired it [25 27, 39]. Other authors have recently reported on RCSJ-based analysis of the three- 
junction Josephson system presenting their own terms and analysis [38]: in these papers, the pumping by an ac drive 
was used in order to generate low frequency modulations of phase and energy and the conclusions confirm a nonlinear 
phenomenon reported previously in simulations and experiments [25-27, 39], namely that it is possible to phase-lock 
Josephson systems by superharmonic pumping. The quantum effects on the other hand should not be sensitive to 
the superharmonic drive and therefore this difference could constitute a discriminant between classical and quantum 
effects. Two of us (MC and NGJ) were co-authors in publications[25, 26] in which the superharmonic pumping of 
phase qubit was observed over a broad temperature range: we recall indeed that nonlinear RCSJ dynamics can model 
a Josephson systcirn over broad parameter ranges and in particular at very low temperatures. Recent observations [34] 
indicate that this is indeed is the case down to lOmK where striking experimental observations of nonlinear effects are 
observed ; therefore it might not be easy to engineer a Josephson system on which all the nonlinear effects are purged 
out. It is also worth noticing that it is not straightforward to draw conclusions from experiments with superharmonic 
pumping of a superconducting circuit due to the varying responses that both the system and the microwave apparatus 
may have at different frequencies. 

The possibility that the behavior of a physical system thought to demonstrate evidence of macroscopic quantum 
superposition can be interpreted with a different approach is not a unique characteristic the Josephson effect. Roughly 
three decades ago scientists were looking for macroscopic quantum effects in charge density-wave systems [40] and it 
was found that the reported results could also be explained by alternative and more "classical" models [41]. It was 
earlier speculated that these kind of arguments would not find space and motivation in Josephson systems, but while 
our findings for resonant switching, Rabi-oscillations, Ramsey-fringes, and spin-echo do not preclude observations of 
quantum behavior in Josephson systems, the evidence of recent years has shown that many observed phenomena in 
this class of systems can be attributed to the RCSJ dynamics developed for Josephson variables in a classical limit [24]. 
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